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ABSTRACT 


Parameter identification is one of the most important 
problemgin groundwater system analysis. For a local problem 
the aquifer parameter can be found by conducting pumping 
tests. But for a regional problem the pumping test method 
will be too expensive, A method using the multicell approach 
is studied to find the aquifier parameters for a region. The 
problem can be stated as an optimization problem, linear or 
non-linear. The ground water flow is represented in bhe 
form of equations, called as balance equations which are 
obtained from the law of continuity and Darcy law. The 
equations are converted into linear constraints to solve 
the optimization problem. Several criteria can be thought 
of in formulating the objective function. The basic principle 
to form the criteria is to minimize the discrepancy occuring 
between two sides of the flow equation. Criteria used to 
formulate the objective function are linear criteria A,B,C,D, 
In criterion E the square of the discripencies is considered. 
Using these criteria an optimal set of values of aquifer 
parameters are found, A method to modify the objective 
function for criterion E is also explained under the name 
of 'Modified Optimization* to get improved results. The idea 
behind the method is to use the first estimates or approximate 
values of parameters so that the results are close to the 
field values. 


X 


The data used for dif Cerent cells do not have equal 
reliability. Hence weights are assigned depending upon bhe 
number of observations in a cell. The method of applycing 
the weights has also been explained. 

All criteria are applied to the aquifer of different 
shapes with different boundary conditions. The results 
obtained are discussed to evaluate the criteria. 



CHAPTER - 1 


INTRODUCTION 


1 ,1 GENERAL: 

In olden times btea societies used perennial streams 
as source of water supply. With time the demand of water 
increased and they searched for other alternatives and 
learned the use of groundwater as source of water supnly. 

New techniques were developed. These include forecasting, 
evaluation, assessment of groundwater etc. All these have 
great importance to society and have wide application in 
field. The groundwater table may be needed within certain 
limits for a number of purposes. Forecasting helps in 
designing and construction of water- re sources systems 
providing measures for regulating the groundwater regime in 
irrigation areas. In agriculture ttne forecasting is very 
much important since the availability of water to the plants 
through their roots from groundwater depends upon depth of water 
table. Also if the water table is too high it may cause 
water logging and loss of crops. Effect of pumping, recharge 
or groundwater flow has to be knovm to forecast the ground- 

I 

water table. 

To keep the restriction on the groundwater table or 
piezometric surface for different purposes, several constraints 
are imposed on certain aspects like pumpage or recharge. The 
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management problem can not be solved untill we know the 
behaviour of the system, 

1 .2 AQUIFER AND ITS PROPERTIES: 

The permeable geologic formations which can store 
and transmit water are termed as aquifers. The formation 
is considered permeable if the rock material con bams cracks , 
fissures and interconnected pores through which bhe water 
can flow. From analytical point of view the properties of 
the aquifer can be divided into microproperties and macro- 
properties, 

1,2,1 Microproperties: 

These properties of the aquifer can be determined 
from a sample taken from it and testing it in laboratories. 
These are porosity, specific yield and permeability. Work 
on porosity has been done by Torzaghi [l943],nuskat [1946,1, one/ 
Luthin [l966j. 

Specific yiol'-l is a measure of water yielding capacity. 
Methods for determining specific yield were outlined by 
l^nzer [l959]. Though specific yield is a microproperty 
it is equal to the coefficient of storage of an unconfined 
quifer, which is a macroproperty. 

Permeability of an aquifer is a measure of its ability 
to transmit liquid through its network of pores. According to 
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Darcy's la-'/ 

V = ki (1 .1) 

where, V = Superficial discharge velocity 

k = pea rme ability 

i = hydraulic gradient. 

Specific velocity = V/n where n is the fiorosity. The 
coefficient of permeability varies more than any other 
engineering property of the aquifer. It can be measured 
by constant head method[For relatively permeable material] 
or falling head test [for low permeability], Fishel [19-42], 
Values of k are different for different directions. 

Generally its value for horizontal direction is more than 
the value for vertical direction, 

1.2,2 Ilacroproperties of Aquifers: 

These are those properties which can be determined 
for the aquifer as a whole from field tests. Laboratory 
studies can yield information on porosity, specific yield, 
and coefficient of permeability of a sample. The permeability 
variation v/ithin the aquifer can be estimated by taking a 
number of samples taken at different points in the aquifer. 
However, quite often the average values obtained are not 
indicative of the triAe condition of the aquifer. To find 
the permeability of the aquifer field tests are conducted 
which are also known as performance tests. With the help of 
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the field tests the aquifer macroproperties: transmissivity 
or coefficient of transmlssibility (T) and storativity or 
coefficient of storage (S) are measured. These coefficients 
are defined as follov/s: 

Transmissivity is defined as rate of flow of v/ater 

■Z 

(mVday) through a vertical strip of aquifer of unit v/idth 
and extending the full saturation height under unit hydraulic 
gradient [Thois,1 935] , Value of T is equal to k multiplied 
by the thickness of the aquifer. 

The storage coefficient S is defined as the volume 
of the water the quifer releases or logjses from storage 
per unit surface area per unit decline of head, In unconfined 
aquifer the storage coefficient is nearly equal to the specific 
yield, 

1,3 DETERMIIIATJON OF AQUIFER PROPERTIES FROM FIELD TESTS: 

The aquifer performance bests may be steady state 
or unsteady- state test. The groundwater flow is represented in 
the form of mathematical model with certain assumptions and 
simplifications applied to the real flow. This model consists 
of mathematical equations which on solving subject to initial 
and boundary conditions give the parameters of the model. If 
the model is very near to the real system the values obtained 
can be treated as real values. The field tests provide data 
on pumping rate and draw-down at one or more observation wells. 



From these data the aquifer parameters T and S can be 
calculated using equilibrium/ nonequilibrium equations. 

Equilibrium equation [Thiem,l906) is based on 
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assumptions . 

(i) Aquifer is homogeneous, isotropic and infinite in 

aereal extent. 

The discharge well penetrates fully, 

T is constant with time and space. 

Pumping is constant for long time and steady state 
is reached, 

Flov^ is laminar 


(ii) 

(iii) 

(iv) 


(v) 


Value of T is given by 
log(rp/r. ) 

T = 0 — J. 

2n (s^-s^ ) 


( 1 , 2 ) 


where, 0 = discharge at pumping well 

radial distances of observation wells from the 
pumping v/ell 

draw-down in observation wells after the steady 
state has reached. 


= 


^ 2 ' ~ 


The non-equilibrium equation was presented by Theis 
[ 1935 ] and modified by Jacob [1950], The assumptions are: 

(i) Well has infinitesimal diameber J 

(ii) water removed from storage is discharged instantaneously 
with decline in head ; 
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(ill) assumptions (i), (ii), (idi) and (v) as in 
equilibrium equation. 

The non-equilibrium equation is a solution of 
second order differential equation subject to the boundary- 
conditions specified above. The equation is \/ritten as: 


Q 

4 ii: T 




V4Tt 


u 


(1.5) 


where, s = drawdo-wn 


t = time after pumping sbarted 


u = 


r^S 

4Tt 


The data from the pumping tests are analyzed using bhe 
analytical model selected^ usually presented graphically 
in the form of type curves. The detailed description of 
aquifer test equation was presented hy Hantush [1964]. 

As stated already a model is an approximate representation 
of the real system, Tt can never represenb bhe system 
exactly. Hence exact results for real situations can never 
be obtained, However a properly conducted pumping coupled 
with careful analysis can lead to good estimates of the 
parameters of the system. 

1 .4 SYSTEM APPROACH: 

The pumping test methods are confined to a small area 
as they give local information only. They can not be applied 



7 


to the whole of the basin. It is too expensive to conduct 
the pumping tests all over the basin. Moreover the assumptions 
implicit in these tests may not be consistent over the extant 
of a basin. For a regional problem where the area is very 
large other approaches are used, A general reviev^ of this 
problem known as inverse problem or some times referred to 
as 'Calibration of the aquifer model' is given by Kisiel and 
Duck stein [1972]. 


Consider a system configuration v^hich has input as 
known information during past, output as the information 
needed in future. 


Input 


> 


Aquifer 

system 


Output 


Detection 


Identification Prediction 


The black box is the system which transform input into output, 
There are in general three types of problems: 

(i) Detection problem 

(ii) Identification or inverse problem 

(iii) Prediction problem. 

Detection involves determining the inputs to the 
system, given hoth the response of the system and. the system 
outputs. 
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Prediction problem involves determining the output 
of the system, given inputs and the system. The problem 
has unique solution. The identification involves determining 
the parameters which govern the response of the system. 

This is an extremely important problem and generally the 
solution obtained is not unique. We have to choose those 
results which on fitting in groundwater model give outpiit 
consistent v,ath the observed outputs, 

1.4 OBJECTIVE OF THE STUDY: 

The objective of the study is to find out the quifer 
parameters by solving the identification problem. The 
problem is solved using different criterion , The results 
obtained by these criteria^ are compared. The effect of 
different boundary conditions on the criterion selected ^is 
also considered in the study, 

1.5 ORGANISATION OF THE STUDY: 

The study made is reported in follov/ing chapters. 

Chapter 2 deals with literature survey. This includes 
theory of groundwater modelling, detailed study of identifi™ 
cation problem and work done in the past. 

Chapter 3 discusses the detailed theory of the 
approach adopted to solve identification problem. 

Chapter 4 deals with the application of the theory 
studied in Chapter 3 to the cases considered. 
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Chapter 5 deals with the results, discussion of 
results and conclusions. 



CHAPl'ER - 2 


LITERATURE REVTE’.-/ 

2.1 GROUNDWATER MODELLING: 

To study the behaviour of the groundwater flow in 
a basin it is not possible to conduct the cause effect tests 
on the prototype. For this modelling of groundwater flow 
is used. This is the process of simulating the groundvrater 
flow for given situation. A model represents the characteri- 
stics of its prototype. With the help of modelling techniques 
the behaviour of the basin can be studied at very low cost. 

The basic principles and types of models have been discussed 
in the following sections, 

2.2 BASIC PRINCIPLES OF GROUNDWATER MODELLING: 

Groundwater movement is governed by some hydrological 
principles. The basic principles are the Darcy 2a.vr and 
continuity equation. Hagen [l839] and Poiseuille [l&^i-6] 
showed that the velocity of the flow in a capillary tube is 
proportional to the hydraulic gradient, Darcy [1856] observed 
that the velocity of flow through porous media (laminar flow) 
is proportional to the first pov^er of the hydraulic gradient. 

This law represents the macroscopic equivalent of the Navier 
Stokes equation for the viscous flow of groundwater. Viscous 
effects are accounted for completely by Darcy law. Mathematical!] 
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it can be represented as 
V = 1^; I 


(a.i) 


where V = velocity of flow in porous media 
k = coefficient of permeability 
I = hydraulic gradienb. 

The continuiby equation is represented as 


(PV^i I 


-(PVy ) H -3—- (PV^ ) -I 


( 2 , 2 ) 


vfhere ,Vy, = velocity of flow in x,y,z directions 

respectively 

P = mass density of fluid 
b = saturated thickness of aquifer 
q = q04j,y,z,t), amount of water added or 
extracted/area of the system per unit 
time 

S =s coefficient of specific storage, 

Q 

From the above two equations the differential equation 
governing the groundwater flow for two dimensional flow in 
quifer having uniform thickness can be obtained [Jacob, 1950 , 
Hantush, ,1 964, De V^iest ,1965] as 

0^h 3^h S 

a a y 2 a z2 <1 a t 


(2,3) 
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where 

S 

- storage coefficient 



= S , b 

s* 


T 

= transmissivity 


h 

= head. 


The groundwater flow can thus be represented in the 
form of a differential equation and solution of thfise equations 
v/ill explain the behaviour of the groundwater flov/. The 
above equation can be used for solving the different types 
of the problems [Jacob, 1950 J Jacob and Lohman,1952 J 
Hantush, 1964J. 

2.3 DIFFERENT TYPES OP OROUNDWATER MODELLING TECHNIQUES: 

The types of models can be catngorlzed as porous 
media models)^* miscellaneous analog models; elecbrical 
analog models based on similarity between Ohm’s law and 
Darcy's law, and digital computer models for numerical 
solution of aquifer flow equations. 

A brief discussion on those models has been given in 
the following sections, 

2.3.1 Porous Media Models: 

(1) Sand Tank Models: The first groundwater model i-zas used 
by P. Porchheimer. He construe tued a sand tank model for 
study of well flow in Graz, Austria in 1898, This model is 
a scale model of the aquifer with scaled down boundaries and 
modified permeability. These are constructed in v/ater tight 
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boxes of various shapes, like rectangular forms, columns, etc. 

Sand tank models have been widely employed for investigating 

( 

a variety of groundwater flow problems, such as well flow, 
seepage, artificial recharge, dispersion, and sea v/ater 
intrusion. 

(2) Transparent Models" A transparent porous media model 
can be constructed by matching the index of refraction of 
the medium v/ith thab of fluid. The moment of the fluid 
can now be traced by means of dye streams. 

2.3.2 Analog Models : 

Flow through porous media obeys laws that governs 
other physical systems including laminar flow of fluids J 
heat, and electricity. Techniques were developed depending 
upon these similarities, ' Non-electrical analog models are 
described in this section. 

(1 ) Viscous Fluid Model or Hele Shaw Model or Parallel 
Plate Models: The first viscous fluid model was 

developed by Hele Shaw [1S97j to demonstrate flow patterns 
around variously shaped boundaries. This model uses the 
analogy between the groundwater flow and the movement of a 
viscous liquid flowing betv/een two closely spaced parallel 
plates , 

(2) Membrane Models: In this mode!! the surface of ground- 

water is represented by stretched thin rubber membrane, A 
central probe, representing a pumping well, deflects the 
membrane. Deflections are analogous to drawdown and can be 



measured by micrometer or by optical techniques, 

(3) Thermal Models'* A thermal analogy can be developed 
for studies of flov/ in homogeneous and isotropic confined 
aquifers. The flov/ of heat in a uniform body of materials 
satisfies the Laplace equation and hence moves on a potential 
flow system in the same manner as groundwater [Hansen, 1 952 ] , 

2,3.3 Electric Analog Model: 

The flow of an electric current is expressed by 
Ohm’s law as, 



where is the current per unit area through a material 

of specific conductivity ^ and dE/dx is the voltage gradient. 
Equation (2,4) satisfies the Laplace equation and when 
compared v/ith Darcy’s equation 

dh 

V = - k -- (2,5) 

dx 

similarity between the two equations is evident. It can 
be seen that velocity V is analogous to electric current 
hydraulic conductivity k to specific conductivity 
and head h to voltage E, Depending upon the similarity 
explained above the following types of the models can be 
created. 

Continuous systems, for which aquifer properties are 
modeled by an electric conductive medium that is continuous 
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in space and ^discre-te systems in which the aqui^^er properties 
are modeled by an assemblage of discrete electric elements 
forming a network. Following are examples of electric 
analog model. 

(1 ) Conductive Liquid Models: Limited to two dimensional 
steady-state situations, the model is formed by an insulated 
tank filled with an electrolyte such as a dilute solution 

of copper sulphate. Boundaries of the tank are scaled to 
represent the aquifer boundaries. To create equipotential 
surfaces the copper electrodes are immersed in the tardc, 

The models have been applied to investigate a variety of 

seepage conditions, drawdov/n near a well field, 

and regional groundwater flows [Debrine, 1970; Todd and Bear, 

1975]. 

(2) Conductive Solid Models: These modelr are best adapted 

for studies of steady flow in confined aquifers. The principle 
is same as that of liquid conductive models, The only 
difference is that Instead of the liquid medium, here the 
solid medium is used. The conducting paper cut to the slope 

of boundary or conducting paint applied to non-conducting solids 
are used for this model [Sherwood, 1963 J Shestakov, 1968], 

(3) Resistance-Capacitance Networks: These networks are 
typically employed to evaluate confined aquifers under non- 
steady, two dimensional flow conditions. Resistance-capacitance 
networks are versatile in that they can readily study a variety 
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of aquifer conditions as well as extensive aquifers requiring 
a large number of nodes. The technique can be extended to 
three dimensional problem. The limitation on application 
of the analogy involves non-linear conditions of varying 
transmissivity in unconfined aquifers and two fluid flov; 
problems [Walton, 1963 , Skibitzke, 1963], 

4. Resistance Networks: It has the same principle as that 

of resistance capacitance networks. It is composed of an 
electronic analyzer coupled to an analog model and the 
analog consists of array of resistors only. Capacitors are 
eliminated so no storage elements are included. These can 
be applied to steady state cases only [Bauwer, 1962J Hubert, 
1963 and Baturic-Rdbcic I 969 .. 

2,3.4 Digital Computer Models: 

With the developments of digital computers it has 
become possible to study any type of aquifer problem by 
means of digital computer model, Kleinecke [l97lj, Pricket L 
[ 1973 ]. The computers are used to solve the mathematical 
models of the aquifers. The aquifer behaviour is represented 
in the form of differential equations. Solution to these 
equation gives the requisite information. 

Two types of the model used on digital computer are 

(i) Finite difference models 

(ii) Finite element models. 
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(i) Finite Difference Methods: In this method, the aquifer 
is divided into a grid and the flow associated with single 
zone of aquifer is analyzed. The flow equation is based 
on equation of continuity. 

Inflov/ - outflow = change of storage 

This relation along with Darcy lav; for the equation of 
motion yields 


6 9 h a 

dx dx ’ dy 


6h 

(T ■— — ) - Q 

dy 



( 2 . 6 ) 


where T = transmissivity 

S = storage coefficients 
h = head 

0 = net flow removed from aquifer. 

The above equation is written in finite difference form 
and solutions on digital computer are obtained [Weber,1979J 
Finder et al., 1970 *, Freeze , 1969 and Prickett, 1971J. 


(ii) Finite Element Mebhods: The f ini be element techniques 
involves solving a differential equation for groundwater flow 
by means of variational calculus. Guymon [1970], Neuman, S.P, 
[l97l], Finder [l977] ^ Pricket [19751. Aquifer is subdivided 
into ’Finite element' of arbitrary shape and sizes. Parameters 
are kept constant for a given element. But they may vary from 
element to element. Solution of the equation is found by 
minimizing a variational function. 
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2*3.5 Hybrid Computer Models: 

Some bimei the finite difference solubion takes very 
long time to complete the iterations. To overcome this 
difficulty a combination of a digital model and a resistance 
network analog is usedy (And the model is called as hybrid 
computer model, [Morris , 1972 j Vemuri ,1969], This approach 
IS advantageous for solving iteration-intensive problems 
such as non-steady flow in unconfined aquifers. 

2.4 IDENTIFICATION PROBLEM: 

Identification is one of the three basic problems 
of system analysis (i.e. Detection, Prediction, Identification) , 
This is the problem in which the parameters of the system are 
predicted with the help of the given inputs and outputs of 
the system. The behaviour of the system is represented in 
the form of partial differential equations. The general 
review of this problem, known as inverse problem or some 
time referred to as ’Calibration' of the model is given by 
Kiesel and Duckstein [1972]. 

The various types of inverse problems are aimed at 
the determination, respectively of the following elements 
(a) parameters i.e. I.P, type 1; (b) Initial conditions, I.P. 
type 2 J (c) Boundary conditions, I.P, type 3} (d) Input, 
type 4; (e) More than one of the above quantities, type 5. 

Most of the reported work in the area of hydrology 
concerns I.P, type 1, The work on other types have been done 
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among others' by : Lattes and. Lions [ 1969 ] on type 2, 

Phillipson [1971 j on type 3; Moench and Kisiel j.1970j on 
problem close to I.P, type 4, Some of the central approaches 
to the I.P. type 5 in the aquifer equation context include 
those of Cole [1973] ,Emsellem and de Marsily [l97l], Frind 
and Pinter [1973], Hairnes et al, l 1968], Nelson [1968J, 

Neuman [1973] » Vemuri and Karplus [1969J» Scarascia and Ponzini 
[ 1972 ) and Stallman [1956]. 

The solution of inverse problem is a powerful method 
for parameter^ ^Identification Solution of this problem tries 
to find out those estimates of the parameters which on putting 
in model reproduce results, closest to the observed aquifer 
response to the real problem. Inverse problem is generally 
used for identifying parameters only, but it can be used to 
get otVier information also, like boundary conditions, initial 
conditions etc. The method to solve the Inverse problem 
can be divided into two groups i.e, indirect and. direct 
[Neuman, 1973]. The former consists of a repetitive solution 
of the direct problem with successive parameter estimates to 
obtain the closest reproduction of the historic water table 
records. The modification in the aquifer parameter values 
at the end of each iteration can be based upon mathematically 
rigorous procedure or subjective reasoning. In the direct 
method of solving inverse problem, the parameters are treated 
as dependent variables and are solved for directly. 
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Cooley and Sinclair [1976] classified the parameter 
calculation methods in three catagories, 

(1 ) Method that involve direct substitution of measured 
(observed) hydraulic head data into the equation 
assumed to apnroximately govern groundwater flow 
followed by calculation of parameters, either directly 
or through minimization of some function of equation 
of residuals* 

(2) The trial and error fitting of measured heads to heads 
calculated with assumed flow equations, 

(3) Optimization methods that minimize the value of a 
selected objective function composed of some measure 
of the difference between observed and computed heads* 

The first method falls under the category of direct 
method of solving inverse problem whereas the second and 
third methods fall under the category of indirect method 
of solving inverse problem, 

2.5 SOLUTION TO THE INVERSE PROBLEM: 

The methods and application of this methods is discussed 
here in brief, 

2,5.1 Methods of Solving: 

As stated above the inverse problem can be solved by 
indirect or direct method. The indirect method of solving 
inverse problem is trial and error method and consists of 
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(1 ) Assuming certain value for the unKno'wn paramebers* 

(2) Solving the direct problem with this assumed value 
i,e, outputs are obtained, 

(3) Comparing the results obtained, in step 2 with the 
actual observation , and if two do not correspond 
within a certain limit then, 

(4) With the help of a suitable (ad^^ustable algorithm) 
changing the values of parameters in step 1 , 

(5) Repeating steps 2 to 4 until satisfactory similarity 
(using a criterion function) between bhe observed 
and the computed value is obtained, 

In direct method the parameters are considered as 
unknowns, no initial estimation are made. The problem is 
solved with the help of optimization techniques (linear 
programming or nonlinear programming) , 

2,5,2 Criterion Functions'. 

Different criteria can be used for optimality 
If b is true parameter vector and J3 is the esbimabed one, 
then the criterion may be , 

(i) The minimization of some function of (p-b). Since 

b is unknown sufficient a priori knowledge is required 
for good results, 

(ii) Minimization of some function or functional of e =u-u 
i.e, the difference between the measured out put 
(including noise) and the model out put. This error 

can be used because o can be made measurable .u=u(b) ,u=u(3) 
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(iii) The minimization of some functional containing 
the measurable process output and the estimates 
of the state vector and the parameter vector . 

^cont3 criterion is v/idely used using least square 
method. The problem reduces to the variational 
problem of finding 3 so what the functional J(,l) 

IS extremized 
J(i3) == u(b) -* u(ii). 

2.5.3 V/ork on Identification Problem: 

Kleinecke [1971] uses linear programming to obtain 
opbimum values of the aquifer parameters. He imposes a 
rectangular grid on the aquifer and discretizes it. The 
aquifer parameters are represented at nodes, surrounded 
by rectangular area. The criterion function preferred by 
Kleincke is to minimize Z where 

Z = S f/iax I I '(2,7) 

i ^ 

where = error at node i for time k, 

Haimes et al, [1968] formulate the problem of 
estimating aquifer parameters with a non-linear criterion 
function and solve it with multilevel optimization. They 
consider a rather specialized case of a cluster of production 
wells in an infinite ao,uifer. The aquifer is divided into 
N wedge shaped homogeneous regions each enclosing a single 
well, A well js supposed to draw water only from the area 
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enclosed in the wedge, and thus the lines delineating the 
regions act as impervious boundaries for the vrell. If 
Total number of wells = N 

Then the values of parameters have to be found such 
that Z is minimiEed. 


Z 


N 

Z E 
1=1 k =1 



- 


(2.3) 


where 



M = number of pressure observation 
^1 k“ calculated pressure head at well i at time k 
h^ observed pressure head at well i,at time k, 
is obtained in close form by the method of images. 


Vemuri Karplus [1969] look unon the problem of 
parameter estimation of an unconfined aquifer as a control 
problem in distributed parameter systems and use a hybrid 


computer to solve it, 

T. ^ Pm/ N . 5 

If [TCx.y) — -• 1 ^ T(x,y) 


9x 


5 h 


] !- 0 = S{x,y) 


Bh 


Bt 


(4.9) 


then the problem is to compute in a region R, the values of 
T and S and boundary of R such that the functional j is 
minimized. 


J 


11-1 

/ / [ h(x,y,t) - h(x,y,t)]^ dRdt 

•ti R 


( 2 , 10 ) 


Method of subjective calibration was used by Gates 
[1974] . The method is effectively a trial and error procedure. 
A suitable numerical scheme is set up to solve the equation 
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S 

h -[- Q - 


T 


ah 

at 


( 2 . 11 ) 


The unknovm parameters and the recharge and discharge are 
now assigned values suboectively. The equation is solved 
to obtain h at grid points. This calculated h is compared 
with the observed h at the same point. Any of the unknown 
parameters are changed subjectively (without any automated 
algorithm) to obtain a better fit between h and h at as many 
as grjd points as possible. This is the simplest method but 
require considerable knowledge of aquifer to start with. 


A definite improvement over the subjective calibration 
described above was made by Lovell et al , [1972J, A consi- 
derable amount of subjective information about the aquifer is 
still required in this method, but this is now quantified 
and used in a logical way. 

Lmsellem and do Marsily [l97l] developed a method 
which they called *An Automatic Solution for the Inverse 
Problem', to determine T for steady state flow in heterogeneous 
isotropic aciuifer. T is determined as an average value over 
an area J the area defining the scale T is determined in the 
process. They start out by assuming the entire area as 
homogeneous and find a value T such that the norms of the 
error (or the selected criterion function) in computed flow 
is minimum. Subsequently the area is subdivided into smaller 
portions and with the same criterion function, the best T value 
is calculated for each subarea. Process is terminated when 
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no further improvement in the criterion function takes 
place , 

Nelson [ 1968 J developed an energy dissipation method 
to evaluate permeability for a heterogeneous, isobropic 
medium. The drawback of method is assumption of zero storage 
coefficient, Chen and Seinfield [l972j discuss the estimation 
of spatially varying parameters in partial differential 
equations. Cannon [196^] solves a different kind of parameter 
identification problem. He considers heat flow in a homogeneous 
isotropic rod of unknown diffusivity, Jones l1963j also gives 
a method to solve for diffusivity when it is a function of 
time. The method depends upon a closed Corm solution of the 
flow equation and can not be applied to a regional scale* Yeh 
and Tause [1 971 J apply the quasi-linearization method as 
developed by Bellman and Kalaba [1965] to a rather idealized 
and local groundwater pumping problem. This method may be 
an alternative to the graphical type-curve method but is 
not applicable to the regional flovi problem, 

A direct method for the identification of parameter of 
non-homogeneous aquifer was developed by Sagar et al. [l97bj. 

The method does not require iterative solution of the aquifer 
equation. The shape of the surface representing observed 
depended variables (i.e, head) is approximated from measured 
samples by means of various interpolation algorithms. Once 
the various derivatives of the dependent variable are approximate 
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The inverse problems reduces locally to the algebraic equations 
of small dimensions. Aquifer conditions of heterogeneity 
and anisotropy are amenable to this method, 

A non-linear optimization method for aquifer parameters 
estimation was developed by Deepak Kasbyap and ^atisb Chandra 
[l98l]. This is a numerical scheme developed to identify 
storage coefficient transmissibility , recharge coefficient 
and orientation of principal permeability directions. The 
scheme is based upon the constrained minimization of sum of 
the squares of the residuces in the Boussinnesq equation. 

Comparison model method to identify the aquifer 
transmissivities in non-homogene ous. Anisotropic aquifer in 
steady state flow was developed by Giansilvio Ponzini and 
Alfredo Lozej [1982J, The method find out the inter block 
transmissivities referred to the sides of the network block. 

The method is based, on the comparison between the real 
gradients and the ones generated by a 'comparison model' 
v/hose initial transmissivity value,. is arbitrarjly chosen and 
constant throughout the surveyed area. The minimum head 
anomaly criterion and the bottleneck criterion enable one to 
select a value or a set of values of initial transmissivity. 

An algorithm wras developed to reach the final solution. 

2.6 MULTICELL APPROACH: 

In the multicell approach, the area is divided into 
cells. The properties of the cell are represented at the 
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centre of the cell (node). The balance equation is written 
for each cell for each time period. The solution of these 
equations will give bhe value of aquifer parameters. 
Optimization techniques (linear programming or quadratic 
programming) are used to reach the set of optimal values. 

The different criterion to reach the optimal value of the 
parameters are given below. 

(i) Minimize the maximum absolute deviation at the nodes, 
(li) Minimize the sum of the maximum absolute deviations 
for all time periods, 

(ill) Minimize the sum of the maximum absolute deviations 
for all cells, 

(iv) Minimize the sum of the absolute deviations for all 
time periods for all cells, 

(v) Minimize the sum of the square of the deviations for 
all cells for all time periods, 

2,7 MODIFIED OPTIMIZATION METHOD: 

Emsellem and de Marslly [ 1969 ] in their 'automatic 

solution to inverse problem' found that sometimes the results 

produced were unsatisfactory or physically meaningless. They 

thought this was due to source/sink (p. .rr, .) function, 

IfO 

However that is not the only reason why unrealistic 
results are obtained Frind and Finder [1973] have elaborated 
more on intrinsical and mathematical features in the solution 
of inverse problem under steady state and found that the system 
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IS highly sensitive to the hydraulic heads. In fact 
hydraulic heads intervene as gradients (difference of heads 
divided by distance) and small variation in head can produce 
large variation in gradient. This is sperially true for 
the cases where the piezometric surface is relatively flat. 

So there are chances of getting eratic results due to * 

(i) Finite difference discretization. 

(ii) Functional representation of hydraulic head, 

(iii) Interpolation in space and time of field, hydraulic 
head data in order bo assign heads to nodes. 

Kleinecke [l971j used linear programming to solve 

inverse problem. He used the criteria of 

1 

Min Z = s i S (2.12) 

i 0 k 

Kisiel and Ducksfein [1972] commenting on Kleinecke 's approach 

stated that the main draw back of the method is that many of 

the T . and S . do not appear in the optimal solution, which 
1 J 

means that these parameters will have to be zero for Z to be 
a minimum. There is no satisfactory explaination for this 
[add Kisiel and Uuckestejn] , although Kleincke speculates that 
with a data base spread over a longer period of time this 
may not happen. V/hatever the length data, however there is 
nothing in the linear programming algorithm that would ensure 
the presence of all parameters in the optimal solution. An 
improvement to Kleineckls approach was done by changing the 



29 


objective function as ^ 

Min Z = E 2 H ( X 2 ^2 ^2J3) 

1 j k 

Although some improvement was there but still unrealistic 
results were present in the solution. Small changes in 
heads produced big changes in values of T and S whereas the 
system was comparatively insensitive to the source sink 
values. 


Agustin Navarro [l982j suggested not to modify 
head values. But modification in objective function was 
made. And sysfcem was the same as used gy Kleincecke. The 
methedology selected v/as based on restricting the resulting 
values of T and S, so that they were not very different from 
the actual field measurements or office estimation of T and S. 
If the T® and S® were the first choice as field parameters, 
the modified function to minimize, was taken as 


k"!"^ T 

F =EEE [X + d f (i - Aa V 

mod ijic '■ 1,0 ^ ij T® . 

1,0 

S. 

+■22 n - ■} 

ij ' Se ^ -I 

1,0 


( 2 , 14 ) 


where d is an arbitrarily chosen parameter at d = 0^ ^mod 
is same as ^ In this case the resulting T. and S v/ill 
produce a set of residuals that are minimum but T. . and S. 


could not be similar lo T® . and S® 

10 1,0 


With d different from 


k-1- 1 /? 

zero, residuals ' will be higher than it is with d -0 

1,0 

and highest value of X. v;ill happen when d is relatively big, 

J 
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and T. and S 


0 0 

i ^ are very similar to T. . and S. 
1>0 IjG lyd 


Acceptability Criterion: 

0 

Between the bounds T. and S, . similar to T. ^ and 

If 3 ^ /i » 3 ^>3 

,k H /2 


from one side (d and X 


minimum from the other 


1,3 ' ' ■■ ' ' ' "" If 3 

(d = 0) there is a range where the practical solution have 

to be found. Acceptable values of T and S. are those 

1,3 1,3 

for v/hich the residuals are minimum. To do this values of 

d is increased from 0^1,2 And the set of values are 

checked. The values vdiich are in practical range for minimum 
d may be selected as true va3.ues. 


Further Improve men bs: 

The objective can further be improved by assigning 
weights to terms in oboective function depending upon the 
area of corresponding cell or numbers of observations wells 
etc. The weights which depend upon the raliabi3.ity of the 
data in the specific cell have been explained in the next 
chapter. 


CHAPTER - 3 


METHODOLOGY 

This chapter deals with the detailed theorey of the 
approach adopted in present work to solve the identification 
problem. The theory of cell approach, with extension to 
multicell model has been discussed in detail. The procedure 
for applying weights to terms in oboective junction has also 
been explained, 

3.1 SINGLE CELL MODEL: 

The entire area of the basin can be represented as 
a single cell. It is a very approximate model and suitable 
for those conditions only v/here the change in the characteristics 
of the basin with respect to space are negligible. 

Consider the basin given in the Fig, 3.1a and 
Fig. 3.1 b, with impermeable boundaries. The river has a 
discharge 0, The entire area is considered as one cell 
with discharge Q and precipitation, artifical recharge etc. 

The following equation can be v/ritten for the above case, 

t [A (N + R - P) - Q] = AS - h^] (3.1) 

The equation is called water balance equation, 

where N = Natural precipitation (volume/unit area/unit time) 

R = artificial recharge (Volume/unit area/unit time) 

P = pumping (volurae/unit area/ unit time) 
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Q = oubflov/ (volume/ time ) 

= average water table elevation at time 
t and t + t respectively. 

If the boundaries are impermeable, the flov/ thab takes 
place through aquifer can be accounted for, by Darcy lav. 
The area is divided into segments, as shown in Fig. 3. 2. 

The volume of water flowing between tv/o segments 
is given by Darcy law 


Q 

v/he re A 

w 

B 

Put in T 
0 


= k . I . A 
= k.I. w . b 

= area of cross section through which the flow 
takes place 

= width of segment, perpendicular to the direction 
of flow 

= saturated thickness of the aquiCer 
= k.b 
= T.I.w 


Total lnflo^^^ for the segments 


N 

S 

i-1 


T w^ 



where N = number of segments. 

Balance equabion for a single cell with previous boundary 


t[ A (W 1 R-P) “ Q ^ WjL J 

= ^ ® t ■ h 1 







. p. 
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Flow 
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Inflow 
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Flow throu9^ ^ y 
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However this model is very approximate since there are 
spatial variations in aquifer characteristics, 

3.2 MJLTICELL APPROACH: 

This approach is adopbed if the area under consideration 
is very large, and there are appreciable variation in aquifer 
characteristics with respect to space. 


Assumptions : 

(i) The flow is in a confined aquifer or in a phreatic 
aquifer where spatial variations of water table are 
small with respect to the thickness of aquifer, 

(ii) Flow IS two dimensional. 

(ill) Aquifer is isotropic but not homogeneous, 

(iv) No flov/ takes place bet\\!een two points within the 
same cell. 

(v) The average properties of the cell are same within 
the cell. 

The continuity equation is written as 


[T 

ax 


8h 


® r n 

— 1_T — “ H- r-p 

0 y 0 y 



where T 
S 
r 
P 


T(x,y) 

SCx,y) 

r(x,y,t; 

p(x,y,t). 


( 3 . 4 ) 
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r and p are replenishment and v/ithdraul from the aquifer/ 
unit horizontal area. 

In multicell approach the entire area is discretized 
into a number of cells having regular areas (rectangle, 
square etc.). Consider an area sho\m in Fig. 3.3. The 
river has barrier boundaries on three sides and constant 
head boundary (river) on one side. The area has been divided 
into four colls. The average heads are shown in figure 
for different cells at a parbicular time. 

In elevation the cross section of the cells is shovm. 

Consider the small gates at the bottom (shovm as opening). 

Gates are shovm to indicate bransmissibility. 

If the gates are openalsuddenly the flov; v/ill take place 
in the directions shovm by arrows depending upon the head 
difference between tv/-o adjacent cells. The different 
characteristics of bhe cells such as T and S are defined 
ab their center (nodes) Fig, 3,4 . Therefore the number of 
paramebers to be found are tvixcc the number oC cells. 

3.3 BALANCE SQUATIOMS: 

The finite difference form of the balance equation (3.4) 
is obtained in the following manner. 

Discretize the whole area into rectangular cells(Fig,3 .4) .Let 

T. . = Transmissivity of cell (i,j) 

^1-1/2 j “ Transmissivity at the side dividing cells at 
(i,j) and (i-1,a) 
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S . 
i»a 

.k+1/2 

■i»a 


Storativity of cell (i,o) 

= Rate of recharge and per unit area of cell (isj) 
at time period 


= Rate of pumping per unit area of cell (1,3) 

' 1 1 3 


"bVi 

at k time period 

Q I = Discharge in direction x through the side at 
^ i-1/2 

i-1/2 to the cell (1,3) 

^k-^V2 _ Head at mid point between the instants k and ki-1, 

1 p J 

Now writing balance eqi\ation for coll 'i,o) (Fig, 3 . 5 a, 3 . 5 b ), 


at 1 [ 0 1 , -0 I ^ a | - Q„i j 

i 1/2 X J _,. W2 y y 


( - J. ./- - 11-1/2 / j+1/2 y j~ 1/2 


By Darcy law 0 = TV/I 

Put the values of Q ,0 in Equation ( 3 . 5 ) 

^ y 


(3.5) 




...i, As 




[ ]/2 < •!- )/2 


(h^ ,..,~ht J 


<^.yy -yj+1 )/2 tay^.^ )/2 


(r1= T p? .) ^x. Ay. = s, .■'tiffi— - h* 1 

1,3 ^ 1 , 3 ^ ^ i'^'^3 1,3 1 i,i i,d'' 


Taking the transmissivities at centres 






(b) 

Fig.3’5 Flow passing through a cell(I,J) 
in a multicell model 
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- ®1,0 -Ar L \.a - ^li = (Ptd - "ho 


Putting Ay_j + .<:)y3+i = ^ ''‘'yj-i-1/2 


= 2Ayj.i/2 


AXi +6X1-1 “ ^^^±-^/^ 


AXi +AXi.,,.| ^ 2^Xj^n/2 


The equation can be -written after dividing through by Ax^/sy^ 
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- b? 


hj 


-hj 


+ [ i-lj. t- ^ -Al J-> J- . j:.lj3 - ■ 

Ay^-i/a 2 ax.^x..^/2 


+ 


.t , t 

J • J -t 

2Ay, ^y, 


0 + 1/2 




2 A 

^1+1/2 

, t lAt 

- n 



At 



, t _ , t 


2AyjAy^,V2 


1,0 


' <Pi,3 - >^1,3 ) (3-6) 


Another scheme is hased on the CrankWicholsen method, where 
the spatial gradients are computed at the mid point of 
the time interval. 

At the beginning, head = h^ _ 

*t J A 

At the end of time interval, head = hT 

1 $ J 

At the mid point of time interval, head = 


Suffix k ' represents the time instant. Then at mid point 

1 ^ ii^'1 

i,0 2 


The equation (3.6) can now be written as 


k^l 

n,f 


-T 


i,d“1 


k-;1 

+ b? ^ T 
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i,0 i-1,0 




1,0 " 1,0 


-i 

1,0 


1^1 ,0 


■^kt1/2 


k-i-1/2 


1,0 


i,d+1 


+ f. 

1,0 


S,. 


i,d 


_ ,.k+V2 


= P 


1,0 


k-i-1/2 

"^i,d 


or 



ff T + 
m 


.pk+V2 



k+1/2 ^ k+V2 
i,0 i,d 


(3.7) 
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a,b,c,d,e,f are known as LHS coerfxcient of balance equatron 

where, _ . , ^ zo 

, . z. ,.^.-1/2 - v,^^+V2 

(3.7a) 


(3,7b) 


k+1/2 

j^kH/2 

- h^"V2 

1,3 ■ 

^i»3 


^^3-1/2 

k+n/2 
’i, 3 

.k+1/2 

:_£r 

2 ^^1 ''^1-1/2 

kH-1 /2 


-h>5^y2 
1. i j. 

1,3 

^^3-V2 

uk+1/2 

- h^"V2 

^k+1/2 

-l-tJ+l. 

2 

'^^i'!l/2 

2Ayj 


.k+1/2 „ v.k-'-l/2 


2 Ax ^ Ax ^_^^2 
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jk-i-1/2 

1,3 


,k+l/2 

'i,d 


k+V2 

i,3 


, k+1/2 _ . k+1/2 


2 A X . Ax. 


i-i-1/2 


,kHl/2 _ j^k-!l/2 
3-' 1/2 


2 

At 


(3.7c) 


(3.7d) 


(3.7e) 


(3.7f) 


The coefficients of the balance equation, a,b,c,d,e,f are 
calculated for each cell and bal:ince equation for each cell 
can be written in the form of finite difference equation for 
each time period. 
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3 . 4 SOLU'l ION riETHOD : 

For brue values of T and S, the equation (5,7) 
must be satisfied and values of its LHS and RHt" will be 
equal. 

If number of cells = { C 
Time steps = N 

The number of equation = NC 

and unknown (T and S) = 2C 

[Kisiel and Duckstein ,1972], 

The system is called, 

Under determined if MC .d 2 C 
Determined if NC = 2 C 

Over determined if NC 2 C 

If the set of 2C equations is solved and the model is a 
true representation of the aquifer, the results obtained 
will be the true values of S and T, But rarely a model is 
a true representation of the practical conditions so the 
values (from model) are not true values. This discripency 
introduces certain noise in each equation* j^ach equation 
becomes an approximation . Since all equations arc approxi- 
mations, our goal will be to find the value which is nearest 
to the true solution i.e. here the problem becomes a problem 
of finding an optjmal set of values which v/ill satisfy the 
equations most closely. This is done with the help of linear 
or quadratic programming. Certain criteria are fixed to 
find a close from solution. 
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As stated above result values and the actual values 
of the model are different . IjHS values of balance equation 
( 3 . 7 ) will not be equal to the RH3 of the balance equation. 
For a better solution the deviation betv/een these t’*ra 
siites of the equation must be minimum. The different 
criteria to get optimal solution are based upon this 
deviation. 

The problem is reduced to LP or MLP problem with 
optimizing the objective function (linear or quadratic) 
subject to given linear constraints. 

3.5 CRITERIA FOR OPTII'HZATION: 

The following five criteria vrill be considered. The 
first four criterions use linear programming ivhereas the 
fifth criterion uses the quadratic programming. 

3.5.1 Criterion At 

This criterion minimizes the maximum absolute deviation 
between two sides of equation (3.7) for all the cells for 
all time intervals. 

Let the value of maxii/tum absolute deviation be X, Then 
the absolute value of any other deviation occurring for any 
time period for any cell will be less than or equal to X, 
Therefore objective function is 

Min F = X 


or 


Max Z 




(3.B) 
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Subject to 


5 k-ll/Z , k^1/2 

2 «m h,3 S -(p 


k^-1/2 

r. . ) I/. X 


This constraints is equivalent to two constraints. 

^ m , ^^'’-V2 ^ ^ ^ k-i-1/2 k-i-1/2 

y K T -h f S , *'X < T) . - r 

in 1, j i,j ^ ^1, j j 


5 k+1/2 

^ '^m ” , j 

m--1 ’ ^ 


a 1, j 


k-i'1/2 kf1/2 

- X < -(p _p ) 

i,j ^ 


T. , S. , X > 0 ¥ i,j,k . 

-L 5 J ? J 

Number of decision variables = 2C + 1 


(3.9: 


Number of constraints 


= 2NC 


where M - number of time periods 
C = number of cells. 


3.5,2 Criterion B: 


This criterion minimizes the sum of the absolute 


values of maximum deviation for all time intervals. For each 
time interval there are C balance equations. Let be 

maximum absolute deviation corresponding to (k-!-1/2)^^ time 
interval. Then all other deviations for that time interval 

i.n H 1 "Ko Icioc* -l-ViQ'in /^v*> ciniiciT 1/2 


Will be less then or equal to X '' 

1 ^ -'^•^1/2 k-i-1/2 . 

I s ^ f, , S, , - (p, , - r. ytl 


ijJ ijO 


k-i-1/2 k+1/2 , 

h- “ r. ) \ ^ 


V 1,3 


Number of such deviations will be equal to the number of time 
periods. The LP problem is formulated as 
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Min F 


k=0 


or Max Z 


M-1 

1 

k=0 


k 1-1/2 
X 


(3.10) 


S-ub^ect to 

5 

2 

m=1 


g T 

OvYl 


m m 


kH-1/2 

i»0 




5 k 1-1/2 

_ £ g T - f. S. 

,] m 1,0 1,0 


kH.1/2 k^1/2 . ^ k;-1/2 

1,0 - 

(3,11) 

kfl/2 ki-V2 ki1/2, 
X “(P-, n -r, ^ 


®i,0 > ° 

Number of decision variables = 2C + N 
Number of constraints = 2NC. 


3. 5.3 Criterion C; 


This criterion minimizes the sum of bhe absolute 
values of the maximum deviations for all cells. Number of 
maximum deviations will be C as number of cells is C. Leb 
the maximum absolube deviation for cell (i,o) Then 

all other deviations at this cell will having lesser or equal 
values to this. 


5 

2 

m=1 


g T 

m 


s. 


1,0 


f , 1 


■^ijO 


k4l/2 

^i,0 


X. 

1,0 


V k 


(3.12) 

The problem now reduces to an LP problem which states, bo 

find S and T, . with, 

1,0 1,0 
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Subject to 


5 

2 

m=', 


^ k+1 /2 /? 

TTi iyj ^ 


ip3 "1,3 Pi,0 


,ki1/2 _ ^ k-!l/2 
i»3 


2 J71 

rfl=1 


„ .k-J-1/2 


i>J 


'i»3 




T O 

, a » i,3’ 1,3 > 


Number of decision variables 
Number of constraints 


^ i»3>k. 

^ C(2+N) 

= 2NC. 


" r 


k-i-1/2) 


ipl 
(3.15) 


In above four criterion the ob 3 ective function was linear 
in next criterion, quadratic form of the objective function 
is used, 

3.5.5 Criterion 3: 


This criterion minimizes the sum of the squares of 
the deviations for all cells for all time periods, 

"“^ipj deviation for cell i ,3 for time period (k'i-1/2). 


Objective function becomes 


Min F = S 4 


i 3 k 


ifJ 


(3.16) 


Subject to 
5 


,k-H/2 ^ 

-L»3 " ^ipj 

(3.17) 


2 Sm^rn ^‘i,3^ "i,3 = PT. -rf 

m=1 


2 g T -1- f. 

^ m 1,3 


k+1/2 _ _ x^‘'V^ s ■n^5’^'V^ _ J'!:^V2 


ipO 


1,3 Pi, 3 


- r 


ip a 



48 


or 


Min 

F 

= Z 

Z 

X 



1 

3 

i>3 

Max 

Z 

^ - Z 

Z 

X , 
1,3 



i 

3 


Subject to 
5 




2 g T I- t. ' S, 

^ p. rn _ t 1 / 2 o 

Sm ^m ^iio ’ 1, o 


1, a 


k+1/2 k+1/2 


- r . ^ 

i- j G 


1,0 • ^^1,0 i,a ^ 

(3.13) 


T, . s ,, X >0 V 1,0 

J-9 d ^ 9 d J J 


Number of decision variables = 3C. 
Number of constraints = 2CN 


3.5,4 Criterion D: 

This criterion minimizes the sum of the absolute 


values of all deviations for all cells for all time intervals, 
kf1/2 

Let X l 3 e the deviation for cell ( 1 , 3 ) for time period (ki'l/2) 

y J 

Then , 


5 k^1/2 _ kv1/2 


2 Sm ^1 -i ^i 'I ' Pi i ^ f- 

m=1 " ^i ,0 


»0 
(3.1A) 


The problem reduces bo the LP problem as 

k+ 1/2 

Min F = Z 2 S X 
i 3 k 


V j,j,k 


or 


k'r1/2 


Max Z =i - ^ £ L x' 

i 0 k 
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k+1/2 


- 1, ■ c" .. 

(3.18) 


T. S > 0 V i,3,k 

I J > J 


Number of decision variables = (2+N) C 
Number of constraints = 2 NC, 

To take into account the negative values of X, a large 
positive number V’ is added to X, such that the sum is 
always positive. Let 


Y X + Y' 

So X = Y - Y' 

The problem can be written as 


Mn F = S j: 2 - Y')^ 

1 j k ^ 

Subject to 

t g T + j^k+V2 3 k-i-1/2 k 1-1/2 

^ °m m ijij 1,0 i,j 2 -fi.-i - r. . 
m=1 ' 1,0 


- y' 




(3.19: 


T q V k'-1/2 

i,0 ' i»0’ ijJ 


V i,0,k. 


The problem can also be solved by putting the value of 
k+1/2 

X. from equation (3.17) in Equation (3.16), Then 


the quadratic problem can be stated as 

k 1-1/2 


Min F = S E E[ 2 
i 0 k i^-l 


e T 
®m m 


*■ ^i i S -(p^+l/2„ k-i-1/2 ,2 

1,0 '^1,0 ^ijO ' 

(3.20) 
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Subject to 


T. 

I 9 3 


> 0 

P 

<>, o 


V 1, jjk . 


3,e REPRESEMTATJOW OF A^UCFER BOUroARIES: 


The boundary condituon of aquifer may be no flov/ 
boundary, constant head boundary or varylng^boimbary » All 
these conditions arc taken care of by the values of bhe 
coefficients of the balance equation. This is done in 
following manners 


Consider a cell ( 1 , 3 ) at the boundary (Figs, 3.6a and 
3 . 6 b) . Th'-' flow to cell ( 1 , 0 ) v/ill be governed by the 

transmissivities of all cells adjoining this cell. To 
incorporabe the flow from boundary consiRer a fictitious 
cell, beyond the boundary. The head in the fictitious cell 
is h;"’!!'' ' with its trans iissivity as Tt . , 


In case the boundary is a barrier boundary there 
will be no flow across the boundary so the head 

V I y J 
2 

odjjjir ao ji, . * The balonce equation of the coll 

^ 7 J 

(i,j) will not have Tt 

1 9 J 

In case the boundary is a pervious boundary, the value 
of is kept equal to the actual head in fictitious 

cell. This may be constant head or varying head with time. 
Nov/ the corresponding coefficient is calculated using these 
head values and the balance equation for cell (i,j) is 
written. In this case one extra transmissivity value' has' t,Q 
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Barrier boundai y 
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i^i j ~ hj j 
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r ig. J*b 


Pervious boundary 

l^f(k-^1/2)_ Q^.^ual head in fictitious cell 
*>1 

(b) 

Representation of aquifer boundaries 
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iDG calculrftod Jor each fie bilious cells. Hence ilie nuiaben 
of parruiv ■ to rs to bo found will become equal 1co twice iiie 
number of colls plus the number of fictitious cells. 

3.7 APPLICITTON OP lALICHTS IN OBJECTI^/E FUNCTION: 

The data used for the problem consists of the 
dimonsj ons of cells, heads and the pumpage or recharge 
(some times known as source/ sink values). The values of 
heads arc obtained from piozornetric wells located in the 
area. 1 1 is evident that tUe cell v/hich has more number 
of piozomotei'’G located in it, will he having more reliable 
data than others v/hich have less number, oi* the data in cells 
far from the piezometric v/ells will be less reliable, as 
compared to those which are comparatively nearer to the 
piezometers. Using this principle certain berms known as 
’reliability index' have been calculated to represent the 
relative reliability of data available in particular cells, 

Those terms v/ore assigned to the cost coefficient in the 
objiectivo function to have more realistic results. 

The procedure to find these indices is as follows: 

Draw the posi tion of the piozornetric wells in the case study 
area. The area represented by each piezometer P is found 
by drawing Thiessen Polygons (Fig* 3* 7), Wow find the total 
area represented by piezometer and the cells contrihuting 
area to this piezometer. The data of the cell having its 
area inside area (partially or fully^ represented by any piezometer 
will be less reliable if area represented by piezometer is large, 




polygons showing area 
ted by piezometers of 

zr 


(rr 



(NJ CO *4 







i.e. tho reliability oi data is inversely proportional 
to area oi polygon, A number N, which is representative of 
a relative reliability oC the data for area represented by 

f'. 

the piezometer is obtained by dividing total area of basin 
by the area of corresponding polygon. Now relative 
reliability of the polygons is knovm as number N. To find 
the reliability of data in a cell the following procedure 
is adopted. 

Let be the piezometers (P^,P 2 ,P^ in Fig. 3.7). 
The Lhosson polygon are dravm. Reliability index R of cell, 




Area of polygon covered by cell 
Total area of the polygon ^ 


Area ABODE 




Area EFGIII ^ ^2 ''' 


0.18) 


Weightage of cell (i,G) can now be calculated as. 
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APPLICATION 

In this chapter the application ol theory presented 
in Chap bo r 3 is discussed. Different shapes of aquifers 
and different boundary conditions are selected for the 
study „ The data used are based on certain realistic situa- 
tions, 

/-I.1 difperp:nt aquifer conditions; 

Roctnngular Aquifer; 

A six coll model for the aquifer it considered, 

(Fig, A.lajbpc) , The boundary conditions chosen are' 

(i) all the sides are barrier boundaries, 

(ii) two sides have varying head and two sides are 
harrier boundaries, 

(ill) two sides have constant head and the other two sides 
nro barrier boundaries. 

4.1 .2 A Stepped Aquifer; 

A fifteen cell model for the aquifer is considered 
with all sides as barrier boundaries(Fig. 4,2). 

A, 2 DATA USED: 

A. 2,1 Time Periods; 

The heads data are assumed to be available for three 
time periods. The duration for these periods are 125,132 and 



Boffier bounciaty 


i 0 ) 


1 


1 ' 


@ 


( 3 ) 


(a) 


© © I © r 




© © 




i 


© 


© 


j 




P«i vie us 
boundary 


(b) 



Constant 
head boundary 


(c) 


Fig.4'1 Six cell model of the aquifer 
with different boundary 
conditions 



4 


Fig.4’2 A stepp(2d aquifer with 
barrier boundaries 
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for bho study ar© shown in Table 4.3a^ 

4,2 » 5 Loft Hand ^ido Coerficlents s 

Tho coofficionts a,b,c,d,c,f of the equation (3 p7) 

are teriTiod LHS coof ficicnts. The definition sketch of 

these coofficlcnts has been shown in Fig, 4. 4 . The values 

of those coofficionts are calculated by putting values of 
k k’l 1 k"^l/2 

» ^L,C! ' ^^i,j i-n equations 3.7(a) to 3.7(f). 

The values of those coo ff tcionts are assigned, to cislls. 


4.2,6 RHS Coo Cf j cionts of the Balance Equation! 


.k-i-1/2 

■i,a 


The RIIS of the balance equation comprises of p^ 

1 y J 

and tViose are calculated by putting the values of 


and rVV^, The net values assigned bo particular 
colls for parbiculab time period . ‘ . 


i|,2,7 Storage Coofficiont and Transmissivity Values! 

To apply tho modified optimization, estimates of T 
and S 1,0, the approximate value of T and S are required in 
dlfCoront colls* For the colls v/here these values are not 
available reasonable values are assumed and shown in 
Table 4,4. 

4,2,0 V/eighbs for Objective Function! 

The weights which have to be assigned bo various colls 
are calculated from fig. (^5?, 3 . 5? Values of reliability indices 
and veightago factors are given in Table 4,b , The weights 



TIME(k) 


TIME(k.l) 



4*4 Coefficients a,b,d.e.f corresponding 
to the celt(i,j)-definition sketch 
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arc applied to bh.e decision variables in the objective 
funcl/ion^ The objoctive Xunction can be written as 


Min P = Z ^ Iw. (X^‘^V2^2 
i j k ^ ^ i > 3 


(4.1) 


4,2,9 Modiiied Objective Function; 


The modiJiod objective function is written from 
equation (2,1/]) as 


F. 


mod 


k H/2 
Z(X. . 
i j k ^ 


H d [ £i:(l- 

ij 



2 S. 

) --r 2 2(1- 

i 3 - 


( 4 . 2 ) 

TVio expected values of T® . and s! are put in the 

3 1,3 

expression and the coefficients of terms involving squares 

and linear berma of S . and T. , are calculated and given 

1,3 ^>3 

in Tabic 4,6, 


U2 I r 
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RESULTS DISCUSSION AND CO iV'SL'J SIGNS 

Pive di Cferent criteria were presented in Chapter 3. 

The results ol applying these criteria to dafferent aquifer 
situations are given in Table 5,1.1 to 5.4.5. 

5.1 OPTIlviTZATION CRITERIA: 

5.1 .1 Cihterion A minimizes the absolute value of the 
maximum deviati.on occuring in the balance equations. The 
criterion will give reasonable results as long as the absolute 
value of maximum deviation is very high compared to the other 
deviatLonsy ds only maximum deviation appears in the objective 
lurction. The linear prop;rarnming problem is not concerned 
with the values of remaining deviations. However, if the 
reraaiiilng deviations have values marginally smaller than the 
maximum deviation the criterion will fail to give reasonable 
results, 

5,1.2 Cri Lcjrion B minimizes the sum of the absolute values 
of maximum deviation occuring for all time periods. This 
criterion v/ill choose the maximum deviation oncuring for each 
time period and will minimize the sum of their absolute values. 
This criterion is stated explictly time wise but not space 
wise in the objective function. The criterion is likely to 
give reasonable results only if the number of time period is 
high. If tlio problem has one time period criterion A and 


B become identical. 
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Criberion G iniiilniizGs th© niaxiinuiii absolute values 
loT all cells, The OJ’ifcerion is explicit spacecvise but not 
tii/ie wise. More tho number of cells, better win be the 
results with this criterion. In the present study the 
number oC cells is greater than number of time periods. So 
the results with this criterion are likely to be nearer to 
real values of the parameters. If the number of cells is 
one the critorJon A and C will become identical in 
op timization, 

5,1, A Criterion D minimizes the sum of the absolute values 
©•r all deviations. The crj terion is explicit time wise as 
wfell as space wise i.o, all are represented in the objective 
f\mcbion for all time periods. 

5.1 ,5 Criterion B minimizes the sum of the squares of 
the residuals ocouring between LHS and RHS of all balance 
equations. In criterion D the absolute values of the 
deviations wore considered. Some times criterion C and D 

do not have foasiblo solutions* For example when the number 
of cells were increased to 15 cells criteria C and D did 
not have feasible results. This problem was again encountered 
in the case of 6 cell model with constant head on two sides 
and barrier boundaries on the other two sides (5.3,1). 

5.2 DISCUSSION OF RESULTS s 

In deciding which criterion is better, the results 
are compared with field values, given in Table A. 4, 



Prom this Table it is seen that for the 6 cell model 
( vanes from 120 iti'/clciy to 170 m^/day and S varies from 
0,06 bo 0,1 /t 7. The results for three different boundary 
conditions were analyzed and discussed below, 

5.2.1 6 Cell Model with Barrier Boundaries on all Sides: 

'The Table 6,1 .1 shows that when criterion A is used, 

three values for T become zero and two values are high namely 

2 

266,75, 304,40 m /day. Criterion B also gives tv^o values of 
T very high , 24^), 60 , 242,16 m^/day. Criterion D which 
gives reasonable T values gives three S values as zero 
compared to two zero S values in the case of A , C and E, 

p 

For criterJan E one value of T is high, 237*50 m /day. 

Thus, crj berion C gives better results than the 
other criteria. 

On comparing the results In Tables 5,1,1 and 5,1,2 
it is observed that criterion C still gives the best results. 

Prom Table 5,1,3, it is seen that the values are 
better at d ~ 0,004, when weights are used. However if 
these results are compared with those obtained from Table 

5.1.1 than are nearer to the approximate values 

given in Table 5.1,5 but S values are better for S^,S^,Sg, 
Thus while no definite conclusion can be drawn as to which 
criterion is best, it appears that botb criteria, criterion 
C and criterion E with modified optimization at d =» 0,004, 
(using weights) are equally good for this case. 
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5,2,2 6 Cell Model Having Varying Head Boundaries on Two 

Sides andBnrrier Boundaries on the Other Two Sides: 

.1 X results obtained in Table 5,2,1 are analyzed 
it Is seen that criterion A gives non zero values to *12 » 

^5 » In which T 2 and (235,00, 201,44 tn^/day 

respectively) are high and (87,5 m^day) is low, 

GritorJon B gives zero values for two cells for both 
parametors S and T, But ^2»'^5 > 253.50 m^day) are 

higher than maximum value of T, Criterion D gives zero 
values of S to 5 colls compared to two zero values of S 
in criterion C, Also the values of Tp (202,25 m^/day) is 
very high. Criterion E assigns tv/o zero values to T, and 
two zero values to S t,e, same as criterion C, But 
( 2 ^ 3.54 m /day) is high. Therefore criterion C gives the 
bast results where T^,T^,T^ are very near to expected 
values and only two S values are zero. 

Comparing the results of Table 5,2,2, for different 
values of d , it is soon that d ~ O .06 gives a set of 
bettor results. Only one S value has zero value and the 
values of T are almost same for different d. 

If the results obtained from Table 5.2.3 are analyzed 
it is seen that if weights are used the criterion E with 
d « G has better results than those v/ith other values of d. 

As for other values of d deviation in S values from their 
expected values is much larger, T is almost same. However, 
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iC this Meb ifj coinpared with the results obtained with 
cril-ciion C ib is scon that overall the criterion C gives 
the host results for this case, 

5,2,3 6 Coll Model with Constant Head Boundaries on. Two Sides 

and Darner Boundan.es on Remaining Sides* 

Tabic 5,3.1 shows that criterion A gives T2 a high 

p 

value (342,34 m'/day), with T.j and being low (69,20 , 33.42 
2 

m /day) , one S value is zero. Criterion C gives high values 
for Tg (310,08 , 291.85 m^/day) and two T and three S 

values aro zero. Criterion B gives only two S values as 
zero and values ot T deviate less from their expected values 
as compared bo criterion C or A, Criterion. D has no feasible 
solubion. Thus bhe criterion B gives better results. 

If bhe results obtained with criterion B are 
compared with those obtained with criterion E it is seen 
:hat criterion E gives only one value of T and two values 
of S as zero as compared to two values of T and two values 
of S as zero, given by criterion B. Though criterion E gives 
a high value to but bhe number of parameters having values 
near to approximate values js more with this criterion. 

So criterion B gives better results as compared to 
criterion B, The modified optimization (Table 5,3.2) has no 
improvement on estimates of parameters. 

If weights are used (Table 5.3.3) 3-1 d — 0.02 T has 
all values non^zero and S has zero values for two cells only. 



69b 


the results obtained with criterion B it is seen that 
criterion B assigns S value as zero for four cells. But 
modified optimization ( d = 1 '' assigns only two values of S 
as zero. For obher values of S, the values are nearer to 
approximate values for modified optimization, T is assigned 
zero values for 9 cell for both the cases. So modified 
optimization with d - 1 is better than criterion B. 

On analyzing results obtained from Table 5, 4,, 3 i,e, 
using weights, it is seen that all S values are non zero at 
d = 1. Four values of T are in permissible limit. If compared 
with results obtained by modified optimization with d = 1 
more values are non zero. 

So, overall, modified optimization v^fith weights at 
d = 1 gives the best results for this case. 
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5,3 CONCLUSIONS AND SUGGESTIONS FOR FURTI-ffiR WORK? 

5,3.1 Conclusions; 

Following conclusions can te drawn from this study: 

(i) The inverse problem can be solved by direct method 
using criteria A, B, C, D or E, The best results 
are generally obtained v/ith criteria B, C or E, 

(ii) The multicell approach can be used even when some 
information is missing, 

(iii) Modified optimization can be used to improve the 
criterion E, This method requires assumption of a 
value for d, d is given arbitrary numbers starting 
from zero. At d = 0 the criterion is identical to 
criterion E, As value of d increases the part of 
the objective function containing d will increase. 

It has been observed that the values of parameters 
obtained have considerable variations with sma^l 
changes in d. For example , (Table 5,1,2 for the case 
of 6 cell model with barrier boundaries), the value 
of d varies from 0,00 to 0.06 but the values obtained 
for parameters have considerable variations, 

(iv) The weights can he applied to improve the results, 
Weights may depend upon the number of observation 
wells or area of cell or importance of cell. For 
instance, the weightage factor applied to the data 
for different cells can be calculated depending upon 



69d 


the area represented by the observation wells. 

If an observation well represent large area, then 
head data will be less reliable, 

5.3.2 Suggestions for Further Work* 

(i) In the present study the head data was assumed 
available for three time periods. This appears 

to be a limitation. Head data should be taken for 
a large number of time periods to get good results 
as it was observed that many values of parameters 
assume lower boundary values (in the present study 
zero values), 

(ii) A second limitation of the present study is that 
no positive lov/er and upper bounds are fixed for 
S and T values. It is desirable in further work 
to assign such bounds based on field tests, 

(iii) The weights used in the present sbudy depend on the 
number of observations in a cell. The method is 
explained in Article 3,7. Othi=*r methods of assigning 
weights should also bo tried in further studies, 
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Table 4,1 • Head Data for 6 Cell Model of the Aquifer 


SI. 

No. 

Cell 

Peizometeric head h. 

1 , 

varying head boundary 

. in 6 cell 

model for 



At the 

At the end 

At the end 

At ,the end 



begining 

of 1 St 

of 2nd 

end of 



of 1st 

Time Period 

2nd Time 

3rd Time 



Time Period 


Period 

Pe nod 



(m) 

(iJi) 

(m) 

(m) 

1 

1 

225.35 

221 .26 

219.87 

222.12 

2 

2 

219.94 

219.91 

218.12 

218.95 

3 

3 

216.33 

218.56 

216.37 

21 5.78 

4 

4 

219.39 

219.08 

218.26 

218.87 

5 

5 

221 .28 

220.87 

220.1 5 

220.98 

6 

■^6 

223.76 

222.66 

222.05 

223.10 

7 

FT1 

223.30 

221 . 01 

219.62 

221 .85 

8 

FTL 

223.26 

222,16 

221 , 55 

222 . 60 

No be 

: For 

6 cell model 

v^rith barrier 

boundaries on all sides 


the 

head data is 

same as in first 6 colls. 



For 

6 cell model 

with constant 

head boundary the 


head data in FT1 

and PT2 are 223.00 m and 

224.00m 


respectively. 
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Table 4,2 : 

Head Data (h^ . ) 

-T-, 2 

for 15 Cell, 

Stepped 




Aquifer 





SI. ' 
No. 

Cell 

(i»D) 

At the begLning 
of 1 st Time 
Period, (m) 

At the Snd 
of 1st Time 
Period(m) 

At the end 
2nd Time 
Period(m') 

At the end 
of 3rd Tim< 
Period (m) 

1 

(1,3) 

218.35 

217.01 


219.09 

217.48 

2 

(1,4) 

223.55 

221 .26 


219.87 

222.12 

5 

(2,2) 

215.24 

214.01 


21 2 . 92 

214.13 

4 

(2,3) 

217.49 

217.05 


217.52 

216.92 

5 

(2,4) 

219.94 

219.91 


218.12 

218.95 

6 

(3,2) 

21 6.89 

215.97 


213.48 

21 5,72 

7 

(3,3) 

216.63 

217.09 


215.95 

216.35 

8 

(3,4) 

216.33 

218.56 


216.37 

215.78 

9 

(4,1) 

217.01 

216.76 


215.79 

216.41 

10 

(4,2j 

218.32 

218.12 


213.04 

218.06 

11 

(4,3) 

219.15 

210.99 


218,15 

218.71 

12 

( 4,0 

219.39 

219.08 


218.26 

218.87 

13 

(5,3) 

221,68 

220,90 


220.35 

221 .06 

14 

(5,4) 

221 .28 

220.87 


220.15 

220.98 

1 5 

(6,4) 

223.76 

222.66 


222.05 

223.10 
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Table 4.3(a): Avera^-re Puinpiag Rate Per Unit 

(m^/da.y/ m^) 

For 3ix Cel l Mo del 


Cell 


During 
1st Time 


Value of 

During 
ITnd Time 



Period 

Per to 

1 

0.0 

0.0 

2 

0.0 

0.0 

3 

1 .7E-3 

3.7842-3 

4 

1 .1932-3 

2.7802-3 

5 

9.3892-5 

1 .7792-4 

6 

8.9002-5 

1 .4452-4 


fie char ge rate/ unit are a 

During 1st Time period = 2, 2208x1 0""^ i 
During ITnd Time period = 7.7574x10’'^ i 
During Illrd Time period = 5.5370x10*"'^ i 


Area 



During 
Illrd Time 

Period 


0.0 
0.0 
1 . BB4S-3 
1.442E-3 
7.453e-5 
4.753E-5 


■7 p 

r / day/m 

i^/day/m^ 

i^/day/m^ 
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Table 4.3(b): 

Average Pumping 
(m^day/m^) 

Rate Ar>ea for 15 Cell Model 

31. Ro. Cell 


Value of ‘ 


1st Time 2nd Time 3rd Tjme 

Period Period Period 




1 . 

(1,3) 

0.0 

0.0 

0.0 

2. 

(1,4) 

0.0 

0.0 

0.0 

3. 

(2,2) 

0.0 

0.0 

0.0 

4, 

(2,3) 

8.4642-4 

2.01 2E-3 

1 ,06E-3 

5. 

(2,4) 

0.0 

0.0 

0.0 

6. 

(3,2) 

1 .746E-4 

2.8430-4 

9.21 7E-5 

7. 

(3,3) 

7.272E-5 

1 ,669E-4 

8.552E-5 

8, 

(3,4) 

1 .700E-3 

3.7a4E-3 

1 . 884 d~3 

9. 

(4,1) 

9.699E-4 

2.21 6E-3 

1.134E-3 

10. 

(4,2) 

9.471 S-5 

2.394E-4 

1 .324E-4 

11 . 

(4,3) 

0.0 

0.0 

0.0 

12, 

<'4.4) 

1 .1930-3 

2.780E-3 

1 .442E-5 

13. 

(5,3) 

0.0 

0,0 

0.0 

14. 

(5,4) 

9.389E-5 

1 . 779E-4 

7.453E-3 

15. 

(6,4) 

8.900S-5 

1 .445E-4 

4,758E-5 



Approximate Values of T and S from Field Tesxs (T:m^/day) 
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Table 4.5; Values of Weightage Factors Assigned to the Cells 
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Table 4*6(a)» Coefficients of T and T for Ifodified Optimization 
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Table 4, 6(b) : Values of Coefficients of for ^fc>ciified Optimization 
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Table 5.1 .1 : Resialxs for 6 cell model, wixh barrier boundaries on all sides 

(T in m^/day) 
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Table *2 : Results of 6 cell model, with barrier boundaries 

on all sides using modified optimization 
(T: m^/day) 


Cell 

d ^ 

0.02 

d=0. 

04 

d=0.06 


T 

^ S 

T 

s’"' 

T 

S 

1 

0,0 

0,0017 

0.00 

0,0005 

0.0 

0.0001 

2 

132,58 

0.0000 

132.44 

0,0000 

132.55 

0,0010 

3 

66,41 

0,0166 

66.21 

0.0034 

66,33 

0.0030 

4 

156,61 

0.0181 

159.58 

0.0000 

159,40 

0.0000 

5 

238.83 

0,0000 

240.38 

0.0000 

240.30 

0 . 0000 

6 

0,00 

0.0049 

0.00 

0.0004 

0.00 

0.0000 



Table 5,1 ,3 ; Results of 5 cell model, v/ith barrier boundaries on all sides using v/eights 
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Table 5*2.1 : Results of 6 cell model, having varying head boundaries on two sides 

p 

and barrier boundaries on reamining sides (T:m /day) 
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Table 5.2,2 : Results of 6 cell model having varying head 

boundaries on two sides and barrier boundaries 
on remaining two sides us mg modified optimization 
method (T: ni /day). 


Cell 

d 

= 0.04 

d=0. 

06 

d=0.08 



T 

s 

T 

S 

T 

S 

1 

0.00 

0,001 9 

0,00 

0.0004 

0,00 

0.000 

2 

95,77 

0.0000 

97.79 

0.0005 

97.59 

0.0001 

3 

76.49 

0,0067 

78.51 

0.0067 

78.49 

0.0067 

4 

106,12 

0.0000 

164.32 

0,0000 

164,23 

0,0000 

5 

246.02 

0.0000 

245.57 

0.0001 

245.43 

0.0000 

6 

0.00 

0.0027 

0,00 

0.0001 

0.00 

0,0000 

FT1 

132.02 


101 .47 


101 ,41 


FT2 

0.00 


0.00 


0.00 




S3 


Table Results for 6 cell model with varying head 

2 

boundaries on two sides using weight (T:m /day) 


Cell 

d=0 


d==0 

.002 

d=0. 

004 

T 

S 

T 

S 

T 

S 

1 

0.00 

0.0089 

0.00 

0.0070 

0,00 

0.0063 

2 

260,90 

0.0000 

263.17 

0.0000 

265.30 

0,0000 

3 

15.69 

0.0437 

14.67 

0.0380 

13.79 

0.0317 

4 

123.60 

0.0983 

124.02 

0.0780 

124.65 

0.0677 

5 

180.56 

0 . 0000 

178.69 

0.0000 

1 77.55 

0.0000 

6 

0,00 

0 . 0248 

0.00 

0.01 50 

0.00 

0.0074 

FT1 

0.00 




0.00 


FT2 

0.00 




0.00 



Table 5,3.1 * Results for 6 cell model with constant head boundary on two sides, and 

- barrier boundaries on remaining sides (T, m /day) 
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Table 5,3.2 : Results JTor 6 celD model v/ith constant head 

boundary on two sides and barrier boundary on 

remaining sides solution using modified 

2 

optimization (T, m /day) 


Cell 


d:=0.01 


d=0.02 

d= 

0.03 


T 

S 

T 

S 

T 

S 

1 

0,00 

0.0021 

0,00 

0.010 

0.00 

0.0004 

2 

77.33 

0.0000 

76.60 

0.000 

76.88 

0.0000 

3 

52,23 

0.0273 

52.09 

0.0116 52.64 

0.0040 

4 

294.27 

0.01 04 

298.76 

0.000 297.96 

0.0000 

5 

330.73 

0.0000 

333.96 

0.000 

334.47 

0.0000 

6 

177.13 

0.0008 

181 .24 

0.001 

180,82 

0.0000 

FT1 

0.00 

0.00 

0.0 


0.0 


FT2 

265.59 

268.70 

272.94 


273 . 62 
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Table : Re cults for 15 cell inode, with barrier 

boundaries on all sides (T: m /day) 


SI; 

Cell 


Criterion 




No. 

( I , J ) 

. a: . 


B 


8 




T 

S 

T 

S 

T 

S 

1 

(1,3) 

0,00 

0.1526 

0.00 

0,02^1 

0.00 

0.0053 

2 

(1,4) 

626.21 

0,1038 

0.00 

0,0000 

n.oo 

0.0114 

3 

(2,2) 

0.00 

0.0000 

0.00 

0.0992 

0,00 

0 . 0000 

k 

(2,3) 

0.00 

1 .6296 

0,00 

0.11 70 

0.00 

o.nooo 

5 

(2.'0 

0.00 

0.5731 

141 ,37 

0 . 0000 

130,75 

0.0000 

6 

(3,2) 

0.00 

0.0424 

0.00 

0.34 A 1 

0.00 

0.0000 

7 

(3,5) 

481 .45 

0.6108 

242.97 

0.0000 

0.00 

0.0000 

8 

(3,4) 

0,00 

0,3754 

0.00 

0.1 777 

1 .39 

0. Oil 54 

9 

( A ,1) 

810.00 

0.0000 

231 .03 

0.4653 

194.82 

0.001 8 

10 

(4,2) 

708.81 

6.7900 

0.00 

0.1 662 

0.00 

0.2121 

11 

(4,5^ 

0,00 

1 .41 30 

0.00 

■ 0.0397 

15.85 

0.0000 

1 2 

(4, d ) 

0,00 

0,00 

144.81 

0,0559 

1 8957^0 

0,0807 

13 

(5,3) 

1 55.96 

1 .3394 

153.83 

0.0835 

0.0 

0.0000 

14 

(5,4) 

103.14 

1 .2800 

176.96 

0.41 80 

21 9.20 

0.0000 

1 5 

(5,4) 

26.43 

0.931^1 

0.00 

0,0000 

0.00 

0,0300 



Table 5.4.2; Resul bs Xot"' 15 cell model having barrier boundary 


on all sides, using, Modified opbimization 
(T: m^/day) 


SI. 

■ Cell ^ 
- 


d = 1 

d=2 


d=-3 


No, 

T 

8 

T 

S ~ 

T 

s”’ 

1 

(1,3) 

0.00 

0 . 0053 

0.00 

0.0053 

0.00 

0.0053 

2 

(1,4) 

0.00 

0.0090 

0.00 

0.0000 

0.00 

0.0000 

3 

(2,2) 

0.00 

0.0028 

0.00 

0.0000 

0.00 

0.0000 

4 

(2,3) 

0.00 

0 . 0000 

0,00^ 

0.0000 

0.00 

0.0000 

5 

(2,4)176.70 

0.0070 

173.37 

0.0001 

1 71 . 69 

0.0000 

6 

(3,2) 

0.00 

0,0003 

0.00 

0.0000 

0,00 

0.0000 

7 

(3,3) 

0.00 

0.0067 

0.00 

0,0001 

0.00 

0.0000 

8 

(3,4) 

1 .38 

0.0454 

1 .39 

0.0454 

1 ,39 

0.0454 

9 

(-^,1 J1 94.14 

0.0818 

1 09 . 57 

0,0818 

1 82 . 91 

0.0818 

10 

(4,2) 

0.00 

0.0361 

0.00 

0.0004 

0,00 

0.0000 

11 

(4,3) 20,6/) 

0,0098 

20.44 

0,0001 

20,24 

0.0000 

12 

(4,4)199.26 

0.0007 

196.13 

0,0807 

192.24 

0,0807 

13 

(5,3) 

0.00 

0.0039 

0.00 

0.0001 

0.00 

0,0000 

14 

(5,4)220,52 

0.0013 

217.63 

0,0000 

214.67 

0,0000 

15 

(6,4) 

0.00 

0.0000 

0.00 

0.0000 

0.00 

0.0000 



Tatle 5.4,3- H-esulxs for 15 cell model v/ith barrier bo‘and.aries on all sides, usins ^reiirn'ts 
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